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ABSTRACT 

We have investigated steady, radial gas outflows (or winds) from galaxies and the de- 
velopment of thermal instability in the hot gas. In order to see the maximum influence of 
_ the instability on the global structure of the galactic outflows, we study inhomogeneous 

\Q ' comoving flows and the non-linear fate of the fluctuations in the flows. We compare 

OP ■ the results with solutions for homogeneous flows. In the case of supersonic flows, the 

Q>^ ■ global structure of inhomogeneous flows is not much different from that of homogeneous 

^^ . flows. However, detailed investigation shows that the average density of inhomogeneous 

f~^ , flows decreases faster than that of homogeneous flows, because local thermal instability 

removes overdense regions in the inhomogeneous flows and reduces the mass flux. We 
Qh| also find that when the gravity of a galaxy is strong, the cold clouds formed from the 

O ' removed gas are distributed in the galactic halo. In the case of subsonic flows, the form 

-*— > ' of inhomogeneous flows is different from that of homogeneous flows near the regions 

d ■ where the flows terminate. The density rise appearing near the regions where the ho- 

mogeneous flows terminate is not seen in the inhomogeneous flows because the local 
K> , thermal instability decreases the mass flux. The cold clouds formed through thermal 

H I instability in the inhomogeneous flow almost all coast to the same maximum radius. 

Subject headings: galaxies: ISM — ISM: jets and outflows — ISM: kinematics and dynam- 
ics 



1. Introduction 

Some observations suggest that gas clouds in and out of galaxies are related to outflows from 
the galaxies. 

In the Galactic halo, there are clouds of neutral hydrogen with a velocity greater than ~ 
80km s~^; they are called high velocity clouds. There has been an argument about the origin. One 
idea is that they result from outflows from the Galaxy; thermal instability in the flows leads to 
cloud formation (Shapiro and Field 1976; Bregman 1980; Li and Ikeuchi 1992). The other is that 



they are made from metal-poor gas falling from the outside of the Galaxy (Oort 1970; Blitz, et al. 
1999). Recently, however, Richter, et al. (1999) found that the iron abundance in a high velocity 
cloud in the Galactic halo is half of the solar value, which supports the former idea. 

Recent observations of Lya absorption system show that the absorbers are associated with 
galaxies (Morris, et al. 1993; Lanzetta, et al. 1995; Stocke, et al. 1995; Bowen, Blades and 
Pettini 1996; Le Brun, Bergeron and Boisse 1996; van Gorkom, et al. 1996; Chen, et al. 1998). 
For example, the spectroscopy and broadband imaging of galaxies toward 3C 273 indicate that 
the Lya absorbers are not distributed at random with respected to the galaxies (Morris, et al. 
1993). Lanzetta, et al. (1995) found an anticorrelation between Lya absorption equivalent width 
and galaxy impact parameter. Moreover, metals have been found in the absorption system, which 
means that the gas is enriched by the products of stellar nucleosynthesis (Ranch 1998). These 
observations suggest that the absorbers are made from the gas ejected from galaxies. In this case, 
thermal instability in the outflows is expected to play an important role in the formation of the 
absorbers. 

Wang (1995a) investigated cooling gas outflows from galaxies. He studied the outflow form, 
the effects of the galaxy potential, the size of outflow regions, and the efficiency of radiative cooling. 
He found that the hot gas in dwarf galaxies can either flow out as galactic winds, or cool radiatively 
to form clouds. In the latter case, the clouds can escape from the galaxies and may become Lya 
absorbers (Wang 1995b). Recently, Nulsen, Barcons and Fabian (1998) also indicated that the 
ejected gas from dwarf galaxies can give rise to a damped Lya absorber. In contrast, Wang (1995a) 
found that massive galaxies like our own tend to confine the gas; the gas released into the halo 
either cools radiatively or results in a galactic corona. 

However, Wang (1995a) considered only spherical homogeneous flows. (Here, homogeneous 
means that density and radius from the galaxy center have a one-to-one correspondence.) In actual 
galaxies, however, there should be density fluctuations that lead to local thermal instabilities. 
Balbus and Soker (1989) examined the linear evolution of the local thermal instabilities in outflows 
and found that the growth is relatively weak. On the other hand, the non-linear evolution is not 
understood well, although complex gas structures observed in galaxies suggest the fluctuations in 
outflows are highly non- linear (e.g. Jogee, Kenney and Smith 1998). Although there are simulations 
done to study the non-linear evolution of instability (e.g. Kritsuk, Bohringer and Muller 1998), they 
cannot resolve fine structure of gas such as shells of supernova remnants and gas blobs composed by 
mass-loss gas from stars; their scales are typically ~ 10 pc and ~ 1 pc, respectively (Mathews 1990). 
If these fluctuations grow through radiative cooling, the overdensities condense and are decoupled 
from outflows. Thus, the substantial quantities of gas may be removed from the outflows, and the 
form of outflows may be changed. As an instance similar to this one. X-ray observations suggest 
that a large amount of gas is decoupled from cooling inflows in the central regions of clusters and 
the cooled gas is deposited in the wide region (Fabian 1994). If it happens to be the same in 
outflows of a galaxy, the clouds born from the decoupled gas may be distributed in a wide range 
of radius of the galaxy. 



However, this consideration may be too simplistic; the motion of overdense regions relative 
to the ambient gas introduces complication. If gravity induces the relative motion, the overdense 
regions mix may with the ambient gas and cannot cool anymore. Hattori, Yoshida and Habe (1995) 
numerically showed that the motion strongly depends on magnetic fields; if the magnetic pressure is 
strong enough, it goes against the gravity and the overdensities can be sustained. Fujita, Fukumoto 
and Okoshi (1997) estimated that the magnetic pressure could sustain the overdensities in actual 
galaxies. 

In this paper, we consider the non-linear evolution of fluctuations in inhomogeneous outflows 
of galaxies and its influence on the global structure of outflows. Since it is difficult to treat the non- 
linear evolution exactly, we first consider the most extreme case. For this purpose, we investigate 
multiphase comoving flows, in which local thermal instabilities are expected to grow most extremely 
because the motion of overdense regions relative to the ambient gas is ignored. We compare the 
results with those of homogeneous flows. 

Our paper is organized as follows. In §2 we summarize our models for both homogeneous and 
inhomogeneous flows. In §3 we give the results of our calculations and describe the differences 
between homogeneous and inhomogeneous flows. In §4 the fate of the clouds formed in outflows is 
discussed. Conclusions are given in §5. 



2. Models 

2.1. Homogeneous Flows 

The equations of mass, momentum, and energy conservation in a steady, symmetric flow are 

1 d{r'^ pu) 



dr 



(1) 



du I dp GM 

^1~ = — 1 2~ ' (2) 

dr p dr r^ 

T^i^Ki^)— 5Am. (3) 

where p, u, p, are respectively the density, velocity, and pressure of the gas, 7 is the adiabatic index, 
G is the gravitational constant, M{r) is the gravitating mass within r, Up is the proton density, A is 
the cooling function, and T is the temperature. These equations can be rewritten as dimensionless 
equations for the Mach number and the dimensionless temperature (equations [2.14] and [2.15] in 
Wang 1995a). 



2.2. Inhomogeneous Flows 

The model adopted here is based on the multiphase cooling flow model of Nulsen (1986). 
Contrary to the model of Nulsen (1986), we consider outflows with relatively large velocity. 

Using the volume distribution function f(p,r,t) at time, t, and position, r, the mass within 
the volume dV, in the density range p to p + dp is given by 

dM = pf{p,r,t)dpdV. (4) 

Differentiating (4), we obtain the equation of mass conservation 

f.V..v/.^.„, ,5, 

where p = dp/dt is the comoving rate of change of the density of gas with density p, and v is the 
velocity. From now on, we assume that the flow is steady and spherically symmetric for simplicity. 

The mass flux function ^ is defined by 

rp 
m{p,r) = / Airr^p'ufdp' , (6) 

where u is the outward radial velocity. As mentioned in §1, we assume that the cooling gas comoves, 
which means that n is a function of r alone. Thus, there is no local relative motion between gas 
blobs before they are decoupled from the flow. From equation (5), it can be shown that ^ satisfies 
the relation 

P^p = -"97 (^^ 

(see Nulsen 1986). Thus, ^ is constant along streamlines in the p, r plane. 

We assume local pressure equilibrium, that is, pressure, p, is also a function of r alone. In 
this case, the streamlines for the cooling gas are determined from two equations; one is the energy 
equation 

pTu'^ = -nlK{T) , (8) 

and the other is the momentum equation 

du _ I dp GM{r) 
dr pdr j-2 ' 

where S is the entropy. The temperature is given by the ratio of p to p (or Up). We assume that 
the cooling function is described by 

A = DT" , (10) 

where D and a are constants. The spatially averaged density is given by 

p=^(cx),r)/ / -^r'^P- (11) 

Jo Pdp 



The total mass flow rate, M, is given by M = Airr'^pu. 
We can integrate equation (8) and give 

{p'-vpf-- = pr^'-\pT^'--^ - p-(2-")) , (12) 

where p is the density at radius r of the gas with density pi at radius ri, p = p(r), pi = p{ri) 
and Pc{r) is the density of the gas at ri which just has infinite density at radius r. The evolution 
of cooling gas along the streamlines is determined by pc- Using equation (12), equation (8) is 
described by Pc 

A T,2 nT" / \ 0.2+0.4Q 



where npi is the number density corresponding to arbitrary density /Oi, and Ti is the temperature 
of gas at pressure pi and density pi . 

The energy and momentum equations can be solved by assuming the functional form of ^. 
According to Nulsen (1986), we consider two cases. One is a power-law form 

10, w > 1 , 

where 

W = (pMn)-^'-") (15) 

for some constant pm- Assuming that A; > 1, we obtain 

P=ip/Plf''jpmil-W,r'/^^~"K (16) 

where 

,^/t + l/(2-a)\ ^^^j 

and 

Wc = (Pc/Pm)-^'-") . (18) 



Thus, equation (9) and (13) can be rewritten as 



d f p\ Jpr, 



Pi J \ ar v^ J 



(19) 



dr \pij pi 

and 

Ul - -c) = (2 - a)^^ U) (1 - ..)--/(2-.) , (20) 

where tcm is the constant pressure cooling time evaluated for r = ri and pi = pm (or Upi = npm if 
represented by number density), that is, 

t - ^^' (21) 
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Since total mass flow rate is M = ^(1 — Wc) , equations (19) and (20) are respectively given by 

n fi / ■ \ -l/fe(2-a) 

d fp\ Jprn fpY'UM] du ^ GM\ .. 



dr \pi J pi \PiJ \ Ml I V dr r^ 

and 
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dr ^ ' tcm \pi) \Mi ^ ^ 



where Mi = M(ri). Using equation (16), we obtain 



M M f p\ I m\ 



(24) 



Attt^P ATrr^Jpra \pij V Ml 

Thus, equations (22) and (23) are simultaneous equations for p and M. 
The other form of ^ is an exponential form, 

^{p,ri)=Ae-'" . (25) 

As is the case of the power-law form (equation [14]), we have the momentum equation 

^(L)^.^(jL-r(j-^9E) (26) 

dr \pij pi \PiJ \ dr r-^ J 

and the energy equation 



where 



^M.-,2-.,l^(ii| , (2, 



'-W^v <^«> 



3. Results 

According to Wang (1995a), we adopt the cooling function for a gas with solar metal abun- 
dance, 

A(T) = 6.2 X 10-i^(T/K)-°-6 cm^ s"! (29) 

(McKee and Cowie 1977; Sutherland and Dopita 1993). This is a good approximation at 10^ < 
T < 4 X 10 K. Note that this cooling function may not be appropriate to study the formation 
of Lya absorbers because their observed abundance is about 1/100 solar. However, although the 
absolute value of the cooling function with the abundance is about one tenth of that of equation 
(29), the temperature dependence is not much different from —0.6 for T < 10^ K (Sutherland and 



Dopita 1993). Thus, the evolution of the flow with solar abundance is similar to that with 1/100 
solar abundance and three or four times larger gas density, because their local cooling times are 
almost identical. 



For a model galaxy, we consider the isothermal mass distribution 

^cir (30) 



GM 2 



with truncated radius rmax- We consider both supersonic and subsonic flows. The Mach number of 
flows, M., is taken as unity at the base of the flows. As discussed by Wang (1995a), it is a natural 
choice at least for the initial condition in supersonic flows, because a smooth transition from a 
subsonic flow at the center to a supersonic flow at large radius requires that M =1 at r = i?, where 
R is the radius beyond which heating ceases (Chevalier and Clegg 1985). The parameters of model 
galaxies are the same as those in Wang (1995a) and are shown in Table 1. Model Dl corresponds 
to a typical dwarf galaxy in which the gas tends to flow out. Model Nl corresponds to a massive 
galaxy like our own in which the gas tends to be confined. Models D2 and N2 are the same as 
models Dl and N2, respectively, but the initial gas temperature and the total mass flow rate are 
larger. For the power-law type of mass flux function (equation [14]), we take k = 2 but the results 
are not much different even if we take k = 10. 



3.1. Supersonic Flows 

As an initial condition, we take Ai= 1 + e (e ~ 10~^ — 10""^) at the inner boundary rj. Since 
the results for models D2 and N2 are qualitatively similar to those for model Dl, we present only 
the results for models Dl and Nl here. In Figures 1 and 2, we show the profiles of pressure, density, 
temperature, velocity, and the ratio of cooling to flow time. For inhomogeneous flows, we plot the 
average density (equation [11]), and the average temperature given by the ratio of the pressure to 
the average density. The cooling time and flow time are defined by 

P 

'^ = (7 - l)n|A(r) (^') 

and 

ti = - (32) 

V 

respectively. In equation (31), we use the averaged density and temperature for inhomogeneous 
flows. In Figures 1 and 2, we also show the total mass flow rate, M, for inhomogeneous flows. 

For homogeneous flows, in model Dl (and D2, N2) hot gas flows out of the galaxy (Figure Id), 
while in model Nl the hot gas is confined, although the gas is first accelerated by pressure (Fig- 
ure 2d). We define a dimensionless variable 

47r(rMir)n2A(T) ' ^''''^ 



That is, 1/x is the fraction of the energy loss due to radiation during the flow time r/vdr for a 
flow with velocity equal to the rotation velocity of the galaxy (Wang 1995a). At the initial radius 
Tj, X ~ 1 in model Nl and x < 1 in others. This means that galactic gravitational potential is 
important in model Nl while it is not important in others, which explains the characteristics of 
the homogeneous outflows. Since initially tf. > t{ in model Nl, the gas first cools adiabatically; 
radiative cooling becomes dominant later on (tc < if)- Note that in model Dl although tc/^f also 
becomes smaller than one for r > 3 kpc for homogeneous flow (Figure le), the gas temperature 
of the flow has already been smaller than 10^ K at r ~ 2 kpc (Figure Ic), that is, equation (29) 
should not be adopted there. 

Broadly speaking, the profiles of pressure, average density, average temperature, and velocity 
of the hot gas for inhomogeneous flows are similar to those for homogeneous flows. However, it is 
shown that the gas pressure and density for the inhomogeneous flows decrease rather faster than 
those for the homogeneous flows (Figures 1 and 2). This is because of the decrease of the mass flux 
(equations [23] and [27]). In fact. Figures If and 2f show that M decreases outwardly, reflecting that 
thermal instability removes overdense regions from the flows. In model Dl, in the region where M 
decreases (r > 1 kpc), the terms of pressure gradient and gravity in equation (9) are small and the 
gas expands freely (u ~ constant). Thus, the decrease of M reduces p because p = M /{Ai^r'^u). In 
model Nl, in the region where M decreases (r > 7 kpc), the gravity term in equation (9) dominates 
the pressure gradient term. This means that the velocity u is determined only by the gravitational 
field and does not depend on the homogeneity of the flow. Thus, the decrease of M reduces p as 
is the case of model Dl. It is to be noted that even in the region where tc/^f > 1, M decreases 
(Figures 1 and 2), because tc is defined for the averaged density and the overdense regions in the 
flows cool ahead of most of the gas. 



3.2. Subsonic Flows 

The initial conditions for subsonic flows are the same as those of the supersonic flows except 
for the initial velocity; we take Ai= 1 — e (e ~ 10~^ — 10"*^) at the inner boundary rj. We present 
the results for models Dl and Nl; the results for other models are qualitatively the same as those 
for model Dl. 

Although gas is confined in a galaxy for all models, the solutions for inhomogeneous flows are 
different from those for homogeneous flows (Figures 3 and 4 ). The density increase appearing 
near the regions where the homogeneous flows terminate does not exist for the inhomogeneous 
flows. Figures 3 and 4 show that tc/if decreases rapidly because the velocity of the flows reduces. 
Thus, radiative cooling becomes important and causes gas condensation near the regions where 
the homogeneous flows terminate. On the other hand, for inhomogeneous flows, radiative cooling 
makes mass-dropout effective and prevents the density increase from appearing. 

In model Nl in §3.1 and in the models discussed in this subsection, radiatively cooled gas 



piles up at finite distances from a galaxy. This appears to be inconsistent with the assumption of 
steady flows. However, if the cold gas turns into dense, compact clouds, the clouds ballistically fall 
back to the galaxy and the gas accumulation does not significantly affect the global structure of 
the outflows unless they are disrupted and mix with ambient gas soon after the formation. The 
stability of the clouds have been studied by many authors. For example, although the clouds are 
subject to Kelvin-Helmholtz instability, Murray, et al. (1993) indicated that the clouds are stable 
for a long time if they are bound by a sufficiently strong (self-)gravitational potential. Moreover, 
Hattori and Umetsu (2000) indicated that the evaporation of the clouds by surrounding hot gas 
may be greatly suppressed by a plasma instability. In the next section, we discuss the motion of 
the clouds including the case that the gravity of a galaxy cannot stop outflows. 



4. Discussion 

We assume that cooled gas in outflows turns to be stable clouds, which move ballistically. The 
initial velocity of the newly formed clouds is expected to be the velocity of the flow at the position 
of their formation if the formation time-scale is small enough. We consider whether the clouds 
escape from a galaxy or fall back under the influence of the galactic potential. 

Since the mass distribution of our model galaxy is truncated at radius rmax; the escape velocity 
is given by 



^Pi'. 



1 + ln 



(34) 



r 

for r < rjiiax- We take rmax = 50 n from now on. If the clouds cannot escape from the galaxy, the 
maximum radius they can reach is given by 



rf = Tc exp 



if A. 

2 \Vax 



(35) 






if rf < rmax, or 

cir 

^^L('^c) - vi 

if rf > rjiiax, where v^ is the velocity of clouds when they form at radius r^ (Wang 1995a). In 
equations (35) and (36), we assume that r^ < rmax- For all the models adopted here, this relation 
is satisfied. 

Figures 5 and 6 show the relation between re and rf for supersonic and subsonic flows, respec- 
tively. First, we discuss supersonic flows. Since in models Dl, D2, and N2, the clouds can escape 
from the galaxy regardless of the position where they are born, we discuss only model Nl. For the 
homogeneous flow, assuming that clouds form when the temperature of gas drops to ~ 10^ K, the 
clouds cannot escape form the galaxy (Figure 5). For the inhomogeneous flows, the clouds also 
cannot escape but rf has a distribution. The clouds that form nearer to the galaxy center have 
smaller rf . By comparing Figure 5 with Figure 2, it is shown that about two third of the clouds 
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formed from the inhomogeneous flows have the orbit which is almost the same as the clouds formed 
from the homogeneous flow (rf ~ 40 kpc); other clouds have smaller apocentric radii. 

Second, we consider subsonic flows. For homogeneous flows, the initial velocity of clouds is 
almost zero, if we assume that they form when the temperature of gas drops to ~ 10^ K. Thus, the 
clouds start to fall back to the galaxy immediately. For inhomogeneous flows, most of the clouds 
cannot escape from the galaxy, although the clouds that form near to n can escape in models Dl 
and D2 (Figure 6). Contrary to the supersonic flows in model Nl (Figure 5), the clouds that form 
nearer to the galaxy have larger rf in models Dl, D2, and N2, because Vc is a strong decreasing 
function of r. However, in model Dl for example, in spite of the wide distribution of rf, it can 
be shown that more than 90% of the gas which originally composed of the outflow can coast only 
to rf ~ 1 kpc by comparing Figure 6a with Figure 3f. This can be said to other models as well; 
most of the gas dropped out of the flow reaches the almost same radius (rf ~ 1 kpc for model D2 
and rf ~ 10 kpc for models Nl and N2). Moreover, the subsonic flows cannot explain the observed 
correlation between galaxies and Lya absorbers because their separations range 12.4 to 157Ah~^ 
(Chen, et al. 1998). 

The clouds formed in the galactic halo fall back toward the galaxy. For our model galaxy, a 
cloud with velocity zero at radius rf has infall velocity. 



■yfaii = -UcirWSln ( ) , (37) 

V Vnall/ 

at radius rfau if rf < rmax- When Vdr = 225 km s~^ and rf = 10 kpc, i;faii = 128 km s~^ at 
Hall = 8.5 kpc. Thus, the observed high velocity clouds may be related to the clouds formed from 
the subsonic flows. 



5. Conclusions 

We study thermal instability in galactic outflows (or winds) and its effect on the global structure 
of the outflows. In order to investigate the maximum effect, we consider comoving inhomogeneous 
(multiphase) outflows in which overdense regions move together with ambient gas; comoving flows 
may be possible if magnetic fields are strong enough. 

We formulate the inhomogeneous flow and apply it to dwarf galaxies and normal galaxies. We 
compared the results with solutions for homogeneous flows. For supersonic flows, the pressure, 
density, temperature, and velocity profiles of hot gas for homogeneous and those for inhomoge- 
neous flows are similar. However, strictly speaking, the density of inhomogeneous flows decreases 
somewhat faster than that of homogeneous flows as the gas expands. This is because local thermal 
instability removes overdensities from the flows owing to local gas condensation. On the other 
hand, for subsonic flows, solutions for inhomogeneous flows are significantly different from those 
of homogeneous flows near the regions where the flows terminate. While radiative cooling and 
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gas condensation raise gas density near the regions where the homogeneous flows terminate, local 
thermal instability in the inhomogeneous flows removes overdensities and prevents an upturn of 
the average density. 

The gas radiatively cooled is most likely to form clouds. The newly formed clouds inherit 
the velocity of a flow at the position of their formation if the time-scale of the formation is much 
smaller than the flow time-scale. For supersonic flows, we flnd that the clouds can escape from 
dwarf galaxies regardless of the homogeneity of flows. On the other hand, for massive galaxies like 
our own, the clouds cannot escape from the galaxies without an extreme starburst; in particular, 
the model of an inhomogeneous supersonic flow predicts that the clouds are widely distributed in 
a galactic halo. For subsonic flows, most of the clouds are conflned in galaxies regardless of the 
homogeneity of flows and the galaxy mass. Most of the clouds formed in an inhomogeneous flow 
reach the almost same radius before being pulled back by gravitation. 

Although the assumption of comoving flows may be too extreme and the effect of thermal 
instability estimated here may be too strong, we think that this possibility cannot be ruled out at 
present. 

I am grateful to K. Okoshi and S. Inutsuka for useful discussions. Comments from an anony- 
mous referee led to significant improvements in the quality of this paper. 
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Table 1. Model Parameters. 
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Fig. 1. — (a) Gas pressure, (b) number density, (c) temperature, (d) velocity, (e) ratio of the 
cooling time and the flow time, and (f ) total mass flux as functions of radius for supersonic flows in 
model Dl. The solid lines indicate homogeneous flows. The dashed and dotted lines respectively 
indicate inhomogeneous flows with the power-law mass-flux function (equation [14]; k = 2) and 
those with the exponential mass-flux function (equation [25]). 
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Fig. 2. — The same as Figure 1 but for supersonic flows in model Nl. 
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Fig. 3. — The same as Figure 1 but for subsonic flows in model Dl. 
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Fig. 4. — The same as Figure 1 but for subsonic flows in model Nl. 
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Fig. 5. — The relation between the position where clouds form and the position where the clouds 
reach for supersonic flows in model Nl. The dashed and dotted lines respectively indicate inho- 
mogeneous flows with the power-law mass- flux function (equation [14]; k = 2) and those with the 
exponential mass-flux function (equation [25]). The diamond shows the result for the homogeneous 
flow. 
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Fig. 6. — The relation between the position where clouds form and the position where the clouds 
reach for supersonic flows. The dashed and dotted lines respectively indicate inhomogeneous flows 
with the power-law mass-flux function (equation [14]; k = 2) and those with the exponential mass- 
flux function (equation [25]). (a) model Dl (b) model D2 (c) model Nl (d) model N2. 



